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ABSTRACT 

Context. Chemical abundances of light elements as beryllium in planet-host stars allow us to study the planet formation scenarios 
and/or investigate possible surface pollution processes. 

Aims. We present here an extension of previous beryllium abundance studies. The complete sample consists of 70 stars hosting planets 
and 30 stars without known planetary companions. The aim of this paper is to further assess the trends found in previous studies with 
less number of objects. This will provide more information on the processes of depletion and mixing of light elements in the interior 
of late type stars, and will provide possible explanations for the abundance differences between stars that host planets and "single" 
stars. 

Methods. Using high resolution UVES spectra, we measure beryllium abundances of 26 stars that host planets and 1 "single" star 
mainly using the A 3131.065 A Ben line, by fitting synthetic spectra to the observational data. We also compile beryllium abundance 
measurements of 44 stars hosting planets and 29 "single" stars from the literature, resulting in a final sample of 100 objects. 
Results. We confirm that the beryllium content is roughly the same in stars hosting planets and in "single" stars at temperatures 
Tcff ^ 5700 K. The sample is still small for T^g < 5500 K, but it seems that the scatter in Be abundances of dwarf stars is slightly 
higher at these cooler temperatures. 

Conclusions. We search for distinctive characteristics of planet hosts through correlations of Be abundance versus Li abundance, age, 
metallicity and oxygen abundance. These could provide some insight in the formation and evolution of planetary systems, but we did 
not find any clear correlation. 

Key words, stars: abundances — stars: fundamental parameters — stars: planetary systems — stars: planetary systems: formation — 
stars: atmospheres 
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1 . Introduction the mixing mechanism a nd rotation behaviour in these stars (e.g. 

^ iPinsonneault et alj|1990h . 

, Since the discovery of the first extrasolar planets, many efforts „ , , r=; 1 r»- A„>,v 7; — , t ^^^. i 

• • , , J : . . • 1 . . J . £ J Several groups (e.g. lGonzalez et al . 2001 : S antos et al.l2001L 

>■ have been made to characterise planet-host stars and to hnd u^^a k i. c j — 7T- Tl — . n , ,,. .' 

, „ , J- • • 1 1 r ^ .1 ^ J . I2004al) have round a correlation between stellar metallicity and 

the features that can distinguish them from stars that do not ■ 1 • 111 

k> , , , " . , — : nPTTrm the presence of giant planetary companion among solar-like 

y\ , have anv known planetarv companion (e.g. lSantos et al.ll2004at , t. u i. 1 1 1 j- i u .u f j - j 

'^«_; r^. — ; ' \, , 1 Urvrvf k rm ^ j r , — '■ — — '^^^ stcus. It has been also largely discussed whether the trend is due 

' IFischer & Valentil I2005D . The s tudy of chemical com p osition ^ ^ • j- v • • • A, f f 1 

and abundances in general (e.g. 'Ecuvillon et al."2004a"S|2Q0i '° Pnmordial ongin, i.e. he frequency of planetary compan- 

. 1 -irvrvr xt . i ^rvrvrv^ j ■ i lon IS a function ot the proto-planetary disc metal composition 

Kjilli et al. 2006; Neves et al. 2009), and, in particular, of the , 17; — / . i ^rvA,^ k • j I . ■ ^ 11 " 

light elements Li. B e and B (e.g. SantosetalJ 120021 [20045 (see e.g.| San tos et al. |2004d), or is due o a posterior pollution , 

rr"^ — 1- . 1 l l ^rvA^I l ^rvrvnk j j ITT ^ . „ I.e. the metalhcity excess IS due to matter accretion after reach- 

llsraelian et al.M2004t 120091) . provided some hints on the innu- ■ ■ / Ti ■ ■ . 1 '^nnn^ »j . . j 

~ p~ — ... J ing the main sequence (see e.g. Pasquini et al. 2007). Most stud- 

ence of a planetary companion in the composition and evo- . „,° • , , , , 

, ^. „ ^ rr, ,. , , J J 1 / N les suggest the primordial scenario although some examples 

lution of stars. These light elements are destroyed by (p,^)- . j. , f, jr 1- ^ 1 lUrvoii i/^Aoo L 



^J^^^^ ^^S^^ ^^^^^'^^^ ^^^^'^y^^J^y (P^^) rnav indicate that both cases occur ( Israelian et al.ll2001l l2()03l; 

reactions at relatively low temperatures (about 2.5, 3.5 and rr—r — nPTTTTT^lT ^ 1 ll -i/-v/-v'? r — 

5.0 X 10^ K for Li, Be and B, respectively). Thus, light ele- l Ashwell et aL Ml Laws et al. ||2003|). 

ments and their abundance ratios are good tracers of stellar inter- However, it is not clear whether a direct relationship be- 
nal structure and allow us to extract valuable information about tween metallicity and probability of planet formation exists or if 

it is only appli cable to the " t ype" o f planets that have been dis 



It IS only appli cable to tne t ype o f planets that have been dis- 
covered so far. ISantos et al.l (l2004ah suggested that two distinct 



* Based on observations obtained with UVES at VLT Kueyen 8.2 m populations of exoplanets can be present depending on whether 
telescope in programme 074.C-0134(A) or not the planet is formed by a metallicity-dependent process. 
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Recently, iHavwoodl (l2009t) has also proposed that the metallic- 
ity excess in the sample of stars hosting planets may come from 
a dynamical effect of galactic nature (i.e. migration of stars in 
the galactic disc), and may not be related with the giant planets 
for mation process . 

llsraehan et al.l (|2004|) found that exoplanet hosts are sig- 
nificantly more Li-depleted than comparison stars but only 
in the range 5600-5850 K. This resu l t has been confirmed 
by IChen & Zhaol (l2006h . iTakeda et alj (l2007h and iGonzalej 
(l2008l) . On the other hand. iGonzaled (12008!) proposed that 
planet-host stars hotter than 5900 K have more Li than the com- 
parison ones. There is observational evidence that Li deple- 
tion is connected with the rotational his tory of the star (e.g. 
iGarcia Loped 11994 iRandich et all Il997h . Thus, similar stars, 
i.e. with similar stellar parameters and age, should have differ- 
ent depletion rates depending on the proto-planetary disc mass 
and composition and its effect on the rotation of the parent star 
( lBouvier2008.) . 

More recently, llsraelian et al.l (l2009l) has confirmed that Li 
is more depleted in Sun-like stars hosting planets than in similar 
stars without detected planets. These stars with and without plan- 
ets, which fall in the effective temperature range 5700-5850 K, 
do not show any differences in the correlation patterns of their 
Li abundances with age and metallicity, pointing to a connec- 
tion between the presence of planets and the low Li abundance 
measured in Sun-like planet-host stars. 

Whatever the metallicity "excess" in planet-host stars is due 
to, light elements can still provide insights on the influence of 
planets on the parent star evolution and in the internal evolution 
of stars in general, they serve as a key for understanding pol- 
lution "events", and also as tracers of the possible planet-disc 
interaction. 

In p revious works directly related with this paper 
dSantos et al. 2002, 2004bx), a comparison between stars host- 
ing planets and stars without known planets revealed that, with 
some exceptions, the two samples follow approximately the 
same behaviour in the variation of Be abundance with the tem- 
perature, similar to Li trend. This behaviour shows a Be abun- 
dance maximum near 6100 K, decreasing towards higher and 
lower temperatures and a Be "gap" for solar-temperature stars. 
ISantos et al.l (12004c') also remarked the possible tendency of the 
planet hosts to be more Be-rich, at least for Teft > 5700 K. 

In this paper we derive beryllium abundances for 26 
stars with planets and one star without planets, and col- 
lect the lithium abundances for most of them. This sam- 
ple enlarge the previous samples of 44 stars with planets 
a nd 29 stars without know n planetary c ompanio ns reported 
in lGarcia Lopez & Perez'de Taorol(ll998h : ISantos et al.l(l2004ch . 
Thus, we end up with a sample of roughly 100 stars where 70 
host planets. We studied the behaviour of these abundances with 
the spectral type and luminosity class in both samples of stars 
in order to further assess the trends found before, and to study 
possible interactions between extrasolar planets and their parent 
stars and its connection to the light-element abundances through 
time. 



2. Observations and data reduction 

We obtained near-UV high-resolution spectra of the targets us- 
ing the UVES spectrograph at the 8.2-m Kueyen VLT (UT2) 
telescope (run ID 074.C-0 134(A)) on 21, 22 December 2004, 
in the blue arm with a wavelengh coverage of /1/13025-3880 A. 
These spectra have a spectral resolution A/ 6 A ~ 70, 000. 
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Fig. 1. HR diagram: surface gravity vs effective temperatures of 
the stars in our sample including the stars in Tables H] |2 ] and |3l 
We overplot solar metallicity tracks from iGirardi et al.l (UOOOl) 
for 0.7 to 1.8 Mq. The luminosity classes of the stars in the sam- 
ple are assigned according with each star position in this diagram 
(See Sect.O. Metallicity is represent by the size of the symbol. 



All the data were reduced using IRAFQ tools in the echelle 
package. Standard background correction, flat-field, and ex- 
traction procedures were used. The wavelength calibration was 
done using a ThAr lamp spectrum taken during the same night. 
Finally, we normalized the spectra by a low-order polynomial fit 
to the observed continuum. 



3. Sample 

The total sample of ~ 100 stars, including 70 planet-host stars, 
cover a wide temperature range from 4500 to 6400 K. Thus, the 
stars hosting planets have spectral types between F6 to K2 and 
three luminosity classes (V, IV and III). The luminosity classes 
were assigned following a criterion based on the position of the 
star in the Hertzsprting-Russell (HR) diagram, and also on their 
surface gravity, i.e. dwarfs (V) have typically \ogg > 4, sub- 
giants (IV) have 3 .3 < log g <4 and giants (III) have log g < 3.3. 
We built an effective t emperature -surface gravity plot with solar 
metaflicity tracks from lGirardi et al. (2000) for 0.7 to 1.8 Mq and 
according to each star position in this diagram (see Fig. [T}, we 
assigned a luminosity class, presented in column 10 of Table [1] 
and column 7 of Tables|2]and [3] In Fig.[T]we represent the metal- 
licity of the objects by increasing the symbol size according to 
three possible metallicity ranges from smaller to bigger size: 
[Fe/H] < -0.5, -0.5 < [Fe/H] < 0.0 and 0.0 < [Fe/H] < 0.5, 
respectively. Thus, our sample of planet-host stars with new Be 
abundance measurements contains 3 giants, 2 subgiants and 2 1 
dwarfs that give a total of 4 giants, 6 subgiants and 60 dwarfs 
in this study. The 30 stars without known planets of our sample 
cover the same spectral types but only two of them are a subgiant 
star. All the new objects are provided in Table \T\ (with 26 stars 
with planets and one star without known planetary companions). 
Stars with and without known planets from literature are listed 
in Tables |2]and [3] respectively. 



' IRAF is distributed by the National Optical Observatory, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under contract with the National Science Foundation. 
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Fig. 2. Best synthetic spectral fits to the observed spectra in the Be ii line region. From upper to lower from left to right, the planet- 
host dwarfs HD 33636 and HD 1 14729, the subgiant HD 99492 and the giant HD 47536 are displayed. 



4. Stellar parameters 

Stellar para meters were taken from the detailed analysis of 
ISousa et aP (12008) when available, otherwise from lSantos et alj 
(|2004a, 2005). These parameters are determined from ioniza- 
tion and excitation equilibrium of Fe i and Fe ii. They used the 
2002 version of the code MOOG0(Sneden 1973) and a grid of 
local thermod ynamical equilibrium (LTE) model atmospheres 
(lKurucdll993 ). The adopted parameters, effective temperature, 
Tgff, surface gravity, \ogg, metallicities, [Fe/H], and masses are 
listed in Table [T] The mean values of the uncertainties on the 
parameters from Santos et al. (2004a, 2005) are of the order of 
44 K for Teff, 0. 1 1 dex for log g, 0.08 km s ' for ^„ 0.06 dex for 
metallicity and the adopted typi cal relat i ve erro r for the masses 
is 0.05 Mq. Whereas from Sous a et aTl (l2008l) these mean un- 
certainties are 25 K for Tetf, 0.04 dex for \ogg, 0.03 km s ' 
for ^, and .02 dex for metallic ity a nd 0.10 for the masses . 
We refer to ISousa et alJ (l2008l) and ISantos et all (l2004ai l2005l) 
for further details. In Tables |2] and |3]we give the main data of 
the sarnp le of stars with and without planets from Sant os et"aD 
( l200 4b c). We note here t he unifo rmity of the adopted stellar pa- 
rameters (see Section 5 in lSousa et al. 2008) . 

One should note that stellar masses are estimated from the- 
oretical isochrones which are strongly sensitive to the adopted 
helium and metal content of the star. In addition, the stellar 
mass, although being a fundamental parameter, is relatively 

The source code of MOOG 2002 can be downloaded at 
|http://verd i.as .utexas .edu/moog .html f 



more uncertain that other stellar parameters, because it depends 
on the effective temperature, surface gravity and metallicity of 
the star. However, unlike the stellar age, the stellar mass for 
main-sequence stars is typically better constrained from theo- 
retical isochrones than fo r subgiant and/or giant stars (see e.g. 
lAllende Prieto et al.ll2004 . 

5. Beryllium and lithium abundances 

The Be abundances were derived by fitting synthetic spectra 
to the data. These synthetic spectra were convolved with a 
Gaussian smoothing profile and a radial-tangential profile to take 
into account the spectral resolution and the macroturbulence, re- 
spectively. The latter was varied between 1.0 and 5.0 km s"', 
between K and F dwarfs (Gray 1992), respectively. A rotational 
profile was also added to account for the projected rotational ve- 
locity, V sin /, of the sample stars, given in Table [1] These were 
calculated from the width of the COR ALIE cross-correlation 
function (see appen dix of iS antos et al. 2002), or taken from 
the literature (mostlv lFischer & Valenti 200 5). A limb-darkening 
coefficient of 0.6 was adopted in all cases, and the overall metal- 
licity was scaled to the iron abundance. 

To derive Be abundances, we fitted all the observed spectral 
range between 3129.5 and 3132 A, where two Be ii lines are lo- 
cated (/13 130.420 and /13 13 1.065 A). We used mainly the Ben 
/13 1 3 1 .065 A line as the feature for deriving the best-fit Be abun- 
dance, and Be ii /13 130.420 A was used only for checking, since 
it is highly blended with other elements lines. 
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Fig. 3. Be abundances as a function of effective tempera- 
ture (top panel) and as a function of mass (bottom panel) 
for the plan et-host stars in ou r samp le plus the sta rs from 
[Garcia Lopez & Perez de Taorol (Il998b : ISantos et al.l (|2002[ 
[2004b c). Filled and open symbols represent stars with and 
without planets. Circles, squares and triangles, are depicted for 
dwarfs (luminosity class V), subgiants (IV) and giants (III). 
In the top panel we superimpose the Be depletion isochrones 
(case A) of Pinsonneault et al. (1990) for solar metallicity and 
an age of 1.7 Gyr. From top to bottom, the lines represent a 
standard model (solid line) and 4 models (dashed lines) with 
different initial angular momentum. We also depict the Be de- 
pletion model (dotted line) including gravity waves provided in 
[Montalban & Schatzman . (.2000.) . 



We adopted the same line list as in ISantos et al.l ( l2004d) . 
which provides a solar photospheric Be abundance of 
logA^(Be) = 1.1 dexB using the same models and tools as in 
this work. The main source of error is probably the placement 
of the continuum. An uncertainty on the continuum position of 
3% yields to an error i n the Be abundance of ~ 0.05 dex for 
solar type stars (see e.g. iRandich et al.ll200"2l) . We estimate this 
uncertainty to be ~ 0.06 dex at Tgff - 5500 K, ~ 0.11 dex 
at Teff = 5200 K, and ~ 0.20 dex at T^b = 4900 K. Another 
weak point is the possible 0.3 dex difference between the solar 
photospheric Be abundance and the meteoritic Be abundance. 
It has been suggested that this could be due to our inability 
to p roperly account for a ll continuum opacity sources in the 
UV dBalachandran & Bell 1998). However, other authors have 
argued that Kurucz atmospheric models are able to reproduce 
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Here we use the notation logMBe) = log[A'(Be)/A'(H)] + 12 
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Fig. 4. Be abundances as a function of effective temperature (top 
panel) and as a function of mass (bottom panel) for only dwarfs 
with (filled circles) and without known planetary companion 
(open circles). Overplotted are the best linear fits found for both 
population together. 



the near-UV absolute continuu m, at least, for stars in the T^g 
range 4000-6000 K ( Allende Pi ieto & Lambert 2000). In fact, 
Balachandran & Bell ( 1998) argue that the Fe i bound-free opac- 
ity should be increased by a factor of 1.6, whi ch is equivalent 
to an increase of 0.2 dex in the Fe abundance. ISmilianic et al.l 
(2009) have tested this possibility for a model of [Fe/H] -0.5 dex 
and they find that the difference in Be abundance is small, 
AlogA^(Be) ~ 0.02 dex. 

For stars with effective temperatures below 5 100 K, the spec- 
tral region surrounding the Be ii /13 1 3 1 A line begins to be 
dominated by the contribution of other species (Mm, C, OH 
lines, etc.) and for the giant stars also the thulium line (Tmii 
/13 13 1.255 A) may play an important role (see lMelo et al.l2005l) . 
This makes it more difficult to fit the observed spectra and thus, 
to measure Be abundance. Although we took into account some 
of these other elements when fitting cooler and giant objects, we 
assume that Be abundances are not so accurate in these cases. 
In addition, the sensitivity of the Be ii to the Be abundance de- 
creases towards lower temperatures. Thus, the total error, includ- 
ing mainly the uncertainties on the effective temperature and 
the continuum location, in the Be abundance is of the order of 
0.1 dex at Teff ~ 6000 K, 0.12 dex at T^g- ~ 5500 K, 0.17 dex at 
Teff ~ 5200 K, and almost 0.3 dex at Teff ~ 4900 K. 

All the Be abundances measured in this work are listed in 
Table [1] We have also added the Be abundance measurements 
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of the 44 sta rs hosting planets and 29 "single" stars from previ- 
ous studies ("Garci a Lopez & Perez de Taoro" 1 9981 ISantos et alJ 
[2002, 2004b c), Hsted in Tables |2] and |3] We have compiled Li 
abundance measur ements, obtained from the Li i /16708 A line, 
from the literature (llsraelian et al.ll2004ll2009b . We also measure 
new Li abundances in four objects (three giants and a subgian t, 
see Table [U in the same way as in (Israelian et al. '2004', '2009J). 
In Fig.|2]we display some synthetic spectral fits to the observed 
spectra of four stars of our sample. 

5. 1 . Ages 

We have gathere d together t h e age s of some of our st ars host- 
ing planets from ISaffe etaP (l2005h . ISaffe et all (l2005h studied 
the correlations between stellar properties with age. They mea- 
sured the chromospheric activity in a sample of 49 stars with 
planetary companions and combining with the literature, they 
obtained age estim ates for 112 obje cts. They applied the cali- 
bration s reported in lDonahue1(ll993b and lRocha-Pinto & Maciell 
(Il998h for the chromospheric activity-age relation but also other 
methods based on isochrones, lithium abundances, metallici- 
ties, and kinematics, and they compared them with the chro- 
mospheric results. They concluded that chromospheric activity 
and isochrone methods give comparable results, and claimed that 
isochrone technique is, in practice, the only tool currently avail- 
able to derive ages for the complete sampl e of planet - host s tars, 
as chromospheric activity calibrations of iDonahue! (1 19931) are 
limited for ages < 2 Gyr and most of the planet-host stars sam- 
ple are > 2 Gyr age. Lithium cannot be used to derived stellar 
ages greater than 1 Gyr Here, we have adopted the ages derived 
from isochrones when available for objects with age < 2 Gyr, 
otherwise we used t hose derived usi ng the chromospheric ac- 
tivity technique with 'Dona hue I (1 19931) calibrations. The adopted 
ages are provided in Table[T] We note that the typical uncertainty 
on the age determination in field stars is very large, between 1 
and 3 Gyr, especially for dwarf stars. 

6. Discussion 

6.1. Be and Li abundance versus effective temperature and 
mass 

In Fig. [3]we plot the derived Be abundances as a function of ef- 
fective temperature (top panel) and as a function of mass (bottom 
panel) for the stars in o ur sample together with those for the stel - 
lar samples r eported in Gar cia Lopez & Perez de Taorol dl998h : 
ISantos et al.l (2002, 2004b cj~ We also distinguish among stars 
with and without planets and among different luminosity classes, 
including dwarfs (V), subgiants (IV) and giants (III). We note 
that we did not find in the literature any determination of the 
stellar mass for all stars in this study, as for instance, no stellar 
mass was available for subgiant and giant planet-host stars with 
Teff < 4900 K. 

Subgiant and especially giant stars have probably changed 
their effective temperatures considerably from their "initial" 
value on the main sequence and their convective envelopes have 
also expanded and reached deeper and hotter regions in the stel- 
lar interiors. Extra mixing may have already occurred and thus, 
the dilution and/or depletion of their light elements do not fol- 
low the normal trend. Therefore, we will take evolved objects 
into account separately (see Sect. l6.1.2] l. 

Considering only dwarfs stars, the addition of the new Be 
abundance of 2 6 planet-host stars in Fig. [3] shows no significant 
changes to what lSantos et alJ (l2004ch found. The overall impres- 
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Fig. 5. Same as Figl3]but for Li a bundance measurements taken 
from llsraeUan et al.f(l2004t l2009h . 



sion is that no clear difference seems to exist between the two 
stellar populations, planet hosts and "single" stars. 

Be abundances decrease from a maximum near Teff = 
6100 K towards higher and smaller temperatures, similarly to 
the behaviour of Li abundances versus effective temperature (see 
Fig. Is). The steep decrease with increasing temperature resem- 
bles the well known Be gap for F stars (e.g. Boesga ard & KingI 
2002), while the decrease of the Be content towards lower tem- 
peratures is smoother, and may show evidence for continuous Be 
burning during the main-sequence evolution of these stars (see 
ISantos et al.ll20045 and references therein). 

In the top panel of Fig. [3] we hav e also posed a set of Yal e 
beryllium-depletion isochrones from iPinsonneault et al.l (Il990l) . 
We depict the depletion isochrones (case A) for solar metal- 
licity and an age of 1.7 Gyr, the standard model and 4 mod- 
els with different initi al angular momentum (see Table 3-6 of 
IPinsonneault et al] ^1990) assuming an initial logA^(Be) = 1.26 
(intermediate value between solar, 1.1 0, and meteoritic, 1.42) 
for all the stars (see lSantos et al.ll2004bl an d references therein). 
As already noticed in Santos et al7('2004b''c). these models agree 
with the observations above roughly 5600 K, but while the 
observed Be abundance decreases towards lower temperatures 
when Teff < 5600 K, these models predict either constant or 
increasing Be abundances. Several possibilities were also dis- 
cussed for this discrepancy between models and Be observa- 
tions: the increasing difficulty to measure Be abundances at low 
temperatures; the lack of near-UV line-opacity in the spectral 
synthesis; the presence of the planetary companions; possible ac- 
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cretion of ^O.SMq by stars in solar neighborhood (iMurrav et al.l 
l200lh : and possible coiTelation with oxygen abundances. 

All of them dismissed as the main possible cause although 
we will further discuss the last one in Sect. 16.41 The better ex- 
planation is probably that these models do not predict correctly 
the behaviour of the Be abundances at the lowest temperatures. 

Mixing by internal gravity waves may provide an 
explanation to t he Li and B e dep letion in cool dwaif 



199U Montalban 1994; 
Charbonnel & Talon. ,2005i) . 



stars ( Garcia L opez & Spruid 
iMontalban & SchatzmanI 119961: 
In Fig. [3] we also depict one model, that includes mixin g 
by gravity waves, from Montalban & Schatzman (■2000). 
ISantos et al.l (|2004b) already pointed out that models includin g 
mixing by internal waves (IMontalban & Schatzman 
reproduce the decrease of Be content in the lower temperature 
regime although still overestimate the Be abundances with 
respect to the observed values. 

In the bottom panel of Fig. [3] we display the Be abundances 
as a function of the stellar masses. The Be abundances show, as 
expected, a similar trend as compared to that of the top panel of 
Fig. [3] However, some stars move with respect to the other stars 
in the figure, specially the evolved stars. As already mention in 
Sect. 16.11 the four giants and one subgiant with Teff < 4900 K 
do not have available mass determinations. A priori, one would 
expect to find the subgiant stars at high masses, but this pic- 
ture changes when having stars with different metallicities. Some 
subgiant stars with planets stay with relatively low masses due to 
their low metal content, whereas other subgiants, like HD 38529 
with a mass of 1.6 Mq and [Fe/H]= 0.4, and with an upper-limit 
Be measurement, have very high masses partially due to their 
high metallicity. In the following, we will concentrate on the Be 
abundances in main-sequence stars. 

6.LL Main-sequence Stars 

In top panel of Fig.|4]we display Be abundances as a function of 
effective temperature only for unevolved stars with and without 
known planetary companion. Here we have excluded the stars 
with Be abu ndances below 0.4 dex which lie in some kind of 
Be-gap (see ISantos et aljr2004bl) . We perform a linear fit, also 
overplotted in Fig. |4] and get a correlation coefficient of 0.74 
with standard deviation of 0.15. If one discriminates between 
the planet hosts and the comparison sample dwarf stars, no sig- 
nificant difference in slope is found when fitting all stars in the 
whole temperature range. 

In Fig. ID (top panel) one can see that whereas in the range 
5700 < Teif < 6200 K, the stars with and without planets mostly 
show similar abundances, at the temperature range 5 100 < T^ff < 
5500 K there seems to be some indication that planet-host stars 
ma y be more Be depleted than "single" stars. 

Ilsraelian et al ](!2009l), based in the Li study of unbiased sam- 
ple of solar-analogue stars with and without detected planets, 
find that planet host have lower Li abundances, being most of 
them upper-limit measurements, than "single" stars, suggesting 
that the presence of planets may increase the amount of mixing 
and deepen the convective zone to such an extent that the Li can 
be burned. 

If this is true, we should see the same effect at lower temper- 
atures in Be burning regime. This may be reflected in the appar- 
ent overabundance of Be in "single" stars respect to planet host 
stars at 5100-5500 K Teff range, in the sense that most of the Be 
measurements in planet hosts are upper-limits. 

However, the sample with Tetf below 5500 K is still too 
small. In addition, as shown in Sect. |5] the expected error on 



the Be abundances in this Teff range is 0.12-0.17 dex. Therefore, 
this calls for new observations to add more Be measurements in 
dwarf stars at these cool temperatures. 

In Fig. |4] (bottom panel) one can see that many points have 
move to a different relative position according to the effective 
temperature due to the different metallicity of the stars with and 
without planets. In particular, the stars without planets which on 
average have lower metallicity appear to have lower masses. So, 
the previous picture with the points well concentrated in a rel- 
atively narrow Teff range, are now spread in a relatively larger 
mass range. In particular, the main-sequence stars with plan- 
ets in the Tf^g range 5100-5500 K have a mean metallicity of 
[Fe/H] ~ 0.18 (with a crpe = 0.22) whereas "single" stars have 
a mean [Fe/H] ~ -0.22 (crpe = 0.20). This explains why the 
slight signature in stars with and without planets seen in the Teff 
plot gets smudged in the log(M/M0) plot. We note that the star 
BD-103166, at T^ff ~ 5320 K and with logA^(Be) < 0.5, is not 
display in Figs.[3]and|4] We did not find any mass determination 
for this star, but we may estimate its mass to be roughly 0.91 Mq, 
and thus having log(M/M0) 0.04. 

On the other hand, the stellar masses trace the stellar posi- 
tions at the beginning of the stellar life and not at their present 
states, which are established by the stellar temperatures. As we 
notice in Sect. 16.31 the stellar metallicity may have an impact on 
the Be burning rate but we do not expect that this explains why 
most of the Be measurements at the Teff range 5100-5500 K 
in planet hosts are only upper-limits. Unfortunately, we cannot 
track the effect of metallicity in the pre-main sequence and main- 
sequence evolution of Be abundances from theoretical models at 
different metallicities since most of the standard models do not 
predict any Be depletion at these cool t emperatures irres pective 
of the metal content of the star (see e.g. lSiess et al ] l2000l) . 

In Fig. |5] we display the Li abundances of the stars in the 
sample versus effective temperature (top panel) and versus mass 
(bottom panel). All the stars in this plot follow the general trend: 
stars with low Be abundances have also their Li severely de- 
pleted. This general trend seen in the top panel of Fig. |5] is less 
clear in the bottom panel due to the larger spread of the Li abun- 
dance measurements in the log(M/Mo) axis. However, there is 
a small group of both stars with and without planets at temper- 
atures Teff < 5670 K which show relatively high values of Li 
abundance, logA^(Li) > 1.2, whereas their Be abundances are 
close to the solar value, logA^(Be) ~ 1.1. Two of these stars are 
dwarf planet-host stars: HD 1237, with a rel atively high chro - 
mospheric activity index, log Rhk - -4.496 dGrav et al]|2006l) . 
classified as an active star, and may be a young object, with 
an age below IGyr; and HD 65216, with a slightly high Li 
abundance, logA^(Li) = 1.28, and a low chromospheiic index, 

logRnx - -4.92. We note that the Sun has logRnx 4.75. 

Among these stars there are also four "single" dwarf stars: 
three of them, HD 36435, HD 43162, and HD 74576, with also 
high chromospheric activity indexes, log Rhk = -4.499, -4.480 
and -4.402, respectively. However, the other object, HD 43834, 
show a relatively low activity (log Rhk - -4.940). This star, 
with a Teff = 5594 K, has a "normal" Be abundance, log A^(Be) = 
0.94 dex and a high Li abundance, logA^(Li) = 2.30 dex, for its 
temperature, according to the general trend. This makes this star 
interestin g in the contex t of pollution. We refer here the studies 
by S antos et all (l2004bllbl) where these four "single" stars have 
been already discussed. 

It is also worth remarking one planet-host dwarf star of our 
sample, HD 142, with a temperature Teif ~ 6400 K, with higher 
Be content than expected from the typical trend, being nearly 
as high as the maximum at T^ff ~ 6100 K. This star also has a 
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cluster with an age of 600 Myr (iBoesgaard & Kinill2002h and in 
the intermediate-age open cluster IC 4651 ('Smi lianic et al .1201 Oh 
with an age of 1.7 Gyr, both with a metallicity of ~ 0.1 dex, very 
similar to the metallicity of HD 142. We note, however, that the 
mass of HD 142, 1.29 Mq, places this star very close to other 
main-sequence stars in the bottom panel of Figs.|3]and|4] 

New observations of Be in stars with and without planets are 
also needed in this temperature range to study possible pollu- 
tion effects. We note that the depletion rates may be different 
for different proto-planetary disc masses and composition, and 
therefore, the rotational history of the star (e.g. Balachandran] 
[T99llaiSretaril200TI) . Thus we will need to find targets with 
similar conditions, i.e. temperature, age, etc. 



logN(Li) 



Fig. 6. Be versus Li abundances of the stars in the sample with 
and without planets. Different colors represent different lumi- 
nosity classes. The stars are splitted into different temperature 
ranges associated with different symbols. 
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Fig. 7. Same as Fig.|6]but for Be abundance versus metallicity. 



higher Li abundance than expected from its effective temperature 
(see Fig.|5]). We also note that the Be abundance maximum may 
need to be slightly shifted towards lower temperatures by one of 
the planet-host stars presented in this work, HD 213240, which 
is located at Teff - 5982 K and have logA^(Be) ~ 1.4. In other 
words, one would have to say that stars with maximum Be abun- 
dances have temperatures in the range 5950 < T^ff < 6150 K. 

The simplest explanation may be that HD 142, classified as 
F7V by S1MBAE0, is just a young star which may be supported 
by its derived and relatively low surface gra vity, log^ - 4.62. 
Howev er, according to the age estimates in Eggenbergeret al. 
(l2007l) . this star has an age of ^2.8-5.93 Gyr which still makes 
lithium content too high. In addition, this star has a relatively 
low c hromospheric activity index log Rhk - -4.853 dGrav et al.l 
l2006h. which may indicate that, in fact, th e age dete r mined 
bv Eggenberger e t al.l (l2007t) may be correct. Grav et ak! (2006) 
classified HD 142 as a chromospherically inactive star. HD 142 
lies in the called Be "gap" according to its derived effective tem- 
perature of Teff ~ 6400 K, but has nearly keep all its Be. In the 
literature one can find stars with similar or even higher Be con- 
tent but typically at T^ff < 6350 K, for instance, in the Hyades 



" The SIMBAD database 
|http://sirnbad.u-strasbg.fr/simbad/| 



6.1.2. Evolved stars 

Following the criterion established in Sect. |3] we have four ob- 
jects clearly classified as giants and eight that could also have 
evolved off the main sequence. 

Five of these subgiants have Li abundance measurements 
and four of them (the planet hosts HD 10697, HD 88133 and 
HD 117176 and the "single" star HD 23249) show relatively 
high values, with respect to the Li trend described by most of the 
dwarf stars. However, these four subgiant stars have Be abun- 
dances consistent with the Be trend shown in dwarf stars (see 
Fig.EJ. 

HD 10697 and HD 117176, as well as the "single" sta r 
HD 23249 were already discussed in Santo s et al.l(l2002[l2004"3) . 



be 



accessed 



at 



These authors discussed the possi bility that pollution by in - 
voking planet engulfment (see e.g. Ilsraelian et al.ll200ll l2003h . 
at least for the two planet hosts, could explain the high Li 
content and the relative "normal" Be content (Siess & Livi"^ 
However, recent models show that many accretion events 
of planetary material could cause Li dep l etion instead of Li 
enhan cement (see e.g. iTheado et al.l 1201 Ol: iBaraffe & Chabrie^ 
I2OIOI) . 

The high Li content and also moderate Be in the metal-rich 
planet host HD 88133, [Fe/H] = 0.33, may be also explained by 
pollution effects. 

The four giant stars, HD 47536, HD 59686, HD 219449, and 
HD 27442, are all of them planet hosts and only present upper 
limits for Li and Be. 



6.2. Beryllium versus lithium 

A beryllium versus lithium diagram can provide information on 
their different depletion rates in main-sequence stars. We have 
already introduced the relationship between Be and Li abun- 
dance in the previous Sect. 16. II but here, we will focus on Fig.|6] 
which depicts the Be abundances versus Li abundances of the 
stars in our sample. In general, there seems to be almost con- 
stant relation between the Be and Li abundances in dwarf stars, 
except for very few dwaifs in the comparison sample. However, 
both stars with and without planets seem to follow the same be- 
haviour. 

In Fig.|6]we also split the Be and Li measurements in several 
stellar temperature ranges. Thus, the objects with Teff > 6000 K 
are situated in the upper right comer, and define the abundance 
maximum, for both Be and Li, and follow a positive correlation. 
The stars with temperatures between 5600 and 6000 K are situ- 
ated in the middle of the diagram, in a nearly flat Be abundance 
through Li abundances ranging from 0.05 to 2.5 dex. The ob- 
jects with temperatures between 5600 and 5100 K are mostly in 
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Fig. 8. Same as Fig.|6]but for Be abundance versus age. 



the left upper-middle panel, also follow an almost constant Be 
abundance at log A^(Be) ~ 0.9 with only upper-limit Li measure- 
ments. This temperature range involve all those stars where Be 
depletion has already taken place and therefore Li is severely 
depleted. Thus, this group shows a Be level lower than that at 
higher temperatures. Nevertheless, 5 objects in this temperature 
regime (two of them are subgiants), both with and without plan- 
etary companions, are located in the right upper side showing 
both high levels of Be and Li abundances. Finally, we separate 
the objects with temperatures below 5100 K as Be abundances 
cannot be determined accurately. These stars can be in turn di- 
vided in three sets: (i) low Be and Li abundances that correspond 
to the giant stars; (ii) low Li abundances but still relatively high 
Be abundances; and (iii) intermediate abundances of both Li and 
Be. The latter contains two objects, one of them being a subgiant 

(see Sect. l6T2b . 

Fig. |6] confirms what found in lSantos et al.l (l2004d) . that Be 
and Li burning seem to follow the same trend for the hottest and 
coolest objects. Li and Be are depleted in a similar way up to 
5600 K, where Be burning stops while Li burning continues up to 
~ 6000 K, where again both Li and Be are destroyed in a similar 
way. Li burning starts to be severe from reff<5900 K while Be 
keeps its low-burning rate until reff~5600 K. Giant stars show 
both Li and Be severely depleted with exception of one object, 
HD 27442, that still keeps a relatively high upper limit to its Be 
abundance. 

In general, Li is more depleted than Be in stars with smaller 
effective temperatures and all stars with and without planetary 
companions show a similar behaviour With the addition of the 
26 planet-host stars and one "single" star, the apparent "lack" 
of planet-host sta rs in the range of Li abundances log A^(Li) ~ 
1.5-2.5, noted bv lSantos et al.l (l2004cl) . has almost disappeared. 
There seems to be roughly the same number of stars with and 
without planets in this range of Li abundances, but this might 
also be due to the small number of stars without planets, 30, in 
this sample, in comparison with the 70 planet-host stars. 

6.3. Be abundance versus /Fe/HJ and age 

In Fig. Q we display the Be abundances versus the metal content 
of the stars with and without planets. We do not see any clear 
correlation but just an increasing dispersion of the Be content 
towards higher metallicities. The differences in the Be content 
are mainly due to the differences in Teff, being lower at lower 
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Fig. 9. Top: Be versus [O/H] abundances in dwarf stars hosting 
planets in the sample and the best linear fit. Bottom: Be abun- 
dance corrected using the linear relation Be-Teff. 



temperatures, but r oughly constant for a ll metallicities. This was 
already noticed by iSantos et al.l (l2004cl) . 

At the highest metallicities, the opacity of the convective 
zone changes considerably, producing a higher rate of Li burning 
and also, although less significant. Be burning. This might ex- 
plain why the Be abundance in planet-host stars has a maximum 
at [Fe/H] ~ 0.20 and then decreases towards higher metallici- 
ties. On the other hand. Galactic Be abundances are know n to in- 
crease with the metallicity (e .g. Rebolo et al., 1988; Molaro et al.l 
Il997t iBoesgaard et al.|[T999b . which in fact is seen in Fig. |7] al- 
though the dispersion in Be abundances and the eiTective tem- 
perature spread is p robably masking this Galactic eff'ect (see also 
ISantos et al.ll2004a, and references therein). 

In Fig. [8] we display the Be abundances versus ages of the 
planet-host stars including those from the literature. We have 
objects in the same temperature regime but with diff'erent ages. 
There is no correlation, indicating that the age, as well as the 
metallicity, is a secondary effect on the level of the Be content 
behind a primary effect driven by the effective temperature of the 
star. However, any interpretation from this plot must be taken 
with caution since the age determinations, specially for main- 
sequence stars, have large error bars and may not be reliable. 

6.4. Be versus [0/H] 

Berylliu m is mainly pro duced in the spallation of CNO nuclei 
(see e.g. iTan et al.]l2009L and references therein). In particular. 
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oxygen gives the largest contribution to the spallation process re- 
sponsible for forming Be. Thus, a comparison between O and Be 
abundances can provide alternative information on the behaviour 
of beryllium. We have collected the oxygen abundances i n stars 
hosting planets of our sample from E cuvillon et al.l (|2006|) . They 
found average oxygen abundances, [O/H], 0.1-0.2 dex higher 
in stars hosting planets with respect to stars without planets. 
However, although planet accretion was not excluded, the possi- 
bility that it is the main source of the observed oxygen enhance- 
ment in planet-host stars is unlikely. Thus, it seems to be related 
to the higher metal content of planet-harbouring stars and both 
Galactic trends of stars with and without planet are nearly indis- 
tinguishable. 

In (top panel) of Fig.|9]we display the beryllium versus oxy- 
gen abundances relative to the Sun only in dwarf planet-host 
stars. We took the [O/H] values derived fr om the near-UV OH 
bands in the spectral range 3167-3255 A (Ecuv illon et al.l2006h . 
because these lines are located close to the Ben 3131 A line 
and therefore would be affected for the same continuum opac- 
ity. We note that the stellar parameters used in this work to de- 
rive Be abundances a re not exactly the same as those used by 
lEcuvillon etan(l2006h to derive oxygen abundances, but in gen- 
eral the differences are very small. We adopted the mea n value 
derive d from the four OH features as given in Ecuvill on et al.l 
(2006). We also note that for some stars in Fig. |3] there is no 
available O measurement. 

The Be content is stratified in different values as the temper- 
ature changes. In the temperature range 5100-5600 K Be abun- 
dance versus [O/H] seems to have a positive slope, just due to 
one dwarf star located at [O/H] ~ -0.25. However, all these 
trends are affected by the strong influence of the effective tem- 
perature on the Be abundances in each star, irrespective of metal- 
licity and oxygen abundance. 

In the top panel of Fig. |9l similarly to Fig. |4l we have also 
displayed the best linear fit, showing a correlation of 0.44 with 
a standard deviation of 0. 18. To study the influence of the effec- 
tive temperature in this Be-O relation, we subtract from each Be 
measurement the Be content given by the Be-Feff linear relation 
in Fig. ID We thus obtain the Be abundance "coiTected" from the 
effect of the effective temperature, logN(Be) - logN{Be{Tfiff)), 
which is shown in the bottom panel of Fig. |9] The Be content 
does not seems to depend much on the initial oxygen abundance 
of the star, as the trend disappears when removing the T^ff effect. 



7. Summary and Conclusions 

We present here new Be abundances of 26 stars harbour- 
ing planets and one star without detected planets, from 
un published UVES spectra, adde d to t he previous s ample 
of 'Garc ia Lopez & Perez de Taorol (Il998l) : ISantos et al.l (l2002l 
[2004b c). We have also compiled Li abundances for these ob- 
jects. The complete sample contains 100 objects, including 
dwarfs, subgiants and giants, with 70 stars hosting planets. We 
study the behaviour of Be abundances of the stars in the sam- 
ple with effective temperature, Li abundance, metallicity, age, 
oxygen abundance and through different spectral types and lu- 
minosity classes. 

In general, the work conclude that the beryllium content in 
stars with and without planetary companion is similar and be- 
haves in the same way. As the burning process is sensitive to the 
mixing mechanism, the presence of a planetary companion is 
not necessary affecting this process in a notorious way, at least, 
at Teff > 5500 K. 



Below Teff ~ 5500 K, the stellar sample is still small but ap- 
pears to be an increase in the dispersion of the Be abundances 
of dwarf stars. In addition, most of the planet-host Be measure- 
ments are not indeed measurements but upper limits. However, 
there is no way to check if the apparent higher scatter of Be con- 
tent at Teff < 5500 K is due to the different stellar metallicities 
since most of the standard theoretical models predict no Be de- 
pletion at these temperatures for any metallicity. 

Be depletion depends on the effective temperature more than 
on the age and metal content. The fact that the stars which host 
planets are richer in different metal contents than "single" stars, 
but have the same Be content, supports the idea of primordial 
origin of these over-abundances in planet hosts. 
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Table 1. Stellar parameters an d Be abundances mea sured in this work. We measured Li abundances for 4 stars. Li abundances 
without a label were taken from llsraelian et aD(l2009b . 
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° Values taken from Sousa et al. 1 200Sj). 
' Values taken from Sant os et al., |20043) . 
° Values taken from Sant os 67^01200^ 

'' Values taken from SIMBAD. 

^ Luminosity class assigned by comparing the stellar parameters with isochrones from lGirardi et allj2000l ) 

^ Values taken from Ecuvillon et al. 1 2006) 

s Values taken from SafFe et al. : 2005) 

' Data taken from Fischer & Valenti 1 200^. 

^ Data taken from Is iaelian et al., ^2004)) . 

^ Data measured in this work. 

* stars without planet companion. 
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Table 2. Stellar paramet ers and Be and Li abundances from literature (iGarcia Lopez & Perez de Taorolll998t ISantos et al.ll2002l 
l2004cl:ISousa et al.ll2008h : Stars hosting planets. 
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° Values taken from lSousa et all 120081) and lSantos et alj)2004ai) 

' Values taken from SIMBAD. 

^ Luminosity class assigned by com paring the stellar parameters with isochrones from lGirardi et al1i2000l) 

Values taken from Ecuvillon et al'f200^ 
' Values taken from Safl'e et al. : 2005) 
^ Values taken from . Garcia Lopez & Perez de Taorol il998l) 
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Table 3. Stellar paramete rs and Be and Li abundances from literature (see lGarcia Lopez & Perez de Taoroll998tlSantos et al.l2002l 
l2004cl:LSousa et al.ll2008h : "Single Stars" 
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HD7615I 


5803 


4.5 


0.14 


1.04 


G2V 


V 


1.02 


1.88 


HD84I17 


6167 


4.35 


-0.03 


1.15 


F8V 


V 


1. 11 


2.64 


HD 186408'' 


5750 


4.2 


0.11 




G1.5Vb 


V 


1.1 


1.24 


HD 189567 


5765 


4.52 


-0.23 


0.87 


G2V 


V 


1.06 


<0.82 


HD192310 


5069 


4.38 


-0.01 


0.80 


K2+V 


V 


<0.6 


<0.2 


HD211415 


5890 


4.51 


-0.17 


0.93 


GOV 


V 


1. 12 


1.92 


HD222335 


5260 


4.45 


-0.16 


0.77 


G9.5V 


V 


0.66 


<0.3I 



Values taken from lSousa et al.||20Ql) and lSantos et alJ|2Q04a|) 

Values taken from SIMBAD. 
^ Luminosity class assigned by comparing the stellar param eters with isochrones from lGirardi et al.] l |2000l) 
^ Values taken from lGarcia Lopez & Perez de Taorol i ll 99 81) 



